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Dark matter halos formed in ACDM cosmologies exhibit a characteristic dependence of 
density on distance from the center. Early studies 2! established pdm ~ or at 
large radii and pdm ~ inside the virial radius. On still smaller scales, the form of Pdm(^) 
was little more than an ansatz since the relevant scales were barely resolved in the A^-body 
simulations. A debate ensued as to whether the profiles were indeed universal, and if so, 
what power of the radius described the dark matter density in the limit r — )■ 0. Subsequent 
studies found both steeper [3|-[6| and shallower [7H9| central profiles. 

The focus of this chapter is the dark matter distribution on sub-parsec scales. At these 
radii, the gravitational force in many galaxies is known to be dominated by the observed 
baryonic components (stellar bulge, nuclear star cluster) and by the supermassive black 
hole. Dark matter densities at these radii are barely constrained observationally; however 
they could plausibly be orders of magnitude higher than the local value at the Solar circle 
(~ IO^^A'/q pc~'^), due both to the special location at the center of the halo, and also to 
interactions between dark matter and baryons during and after formation of the galaxy. 
High dark matter densities make the centers of galaxies preferred targets for indirect de- 
tection studies, in which secondary particles and photons from the annihilation or decay of 
supersymmetric dark matter particles are detected on the Earth. 



I. PHENOMENOLOGY OF GALACTIC NUCLEI 



The distribution of baryonic matter at the centers of galaxies is relevant to the dark 
matter problem for several reasons. 

• Many dynamical processes affect the dark and luminous components in similar ways. 
The distribution of stars at the center of a galaxy can tell us something about the 
distribution of dark matter. 

• If the nuclear relaxation time (eq. |5]) is shorter than the age of the universe, stars will 
exchange kinetic energy with dark matter particles causing the dark matter distribu- 
tion to evolve in predictable ways. 

• Supermassive black holes appear to be ubiquitous components of galactic nuclei. De- 
pending on its mode of growth, a SMBH can greatly increase, or decrease, the density 
of dark matter in its vicinity. 
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Galactic nuclei are the innermost regions of stellar spheroids: either elliptical galaxies, 
or the bulges of spiral galaxies. Most galaxies are too distant for individual stars to be 
resolved, and descriptions of their structure are generally based on their luminosity profiles, 
the surface brightness of starlight as a function of distance from the center. Luminosity 
profiles of galactic spheroids are well fit at most radii by Sersic's [10] law, 

lnJ = lnJo-6i?^/" (1) 

where / is the surface brightness at projected radius R and n is the Sersic index; = 4 is the 
de Vaucouleurs law Sersic's law predicts a space density that increases as p^, ~ ^-in-i)/n 
toward the center, or p^, ~ 7--3/4 ^^^q case of de Vaucouleurs' law. However in the best- 
resolved galaxies - which include both the nearest, and the intrinsically largest, galaxies - 
deviations from Sersic's law often appear near the center. Bright spheroids exhibit deficits 
with respect to Sersic's law, or cores. Faint spheroids exhibit excesses, or nuclear star clusters 
(NSCs), with sizes in the range 1 — 100 pc and luminosities in the range 10^ — IO^Lq 12|, 13|. 



The transition from deficit to excess occurs at spheroid luminosities of roughly IO^^Lq 
While NSCs are generally unresolved (a notable exception [l^ being the NSC at the center 
of the Milky Way), cores in luminous elliptical galaxies can extend hundreds of parsecs. 

Galactic nuclei also contain supermassive black holes (SMBHs). In a handful of galaxies, 
the presence of the SMBH is indicated by a clear Keplerian rise of stellar or gas velocities 
inside a radius ~ rh, the gravitational influence radius: 

Th = GM,/a\ (2) 

Here M, is the black hole mass and a is the one-dimensional, rms velocity of stars in the 
spheroid. In other galaxies, indications are seen of a central rise in velocity but the implied 



SMBH mass is very uncertain 16|. Among the ~ dozen galaxies with well-determined SMBH 



17|: 



masses, there is a remarkably tight correlation between M, and cr, the M, — a relation 

M./lO^Mo ^ 1.66((t/200 km s"^)", a ^ 4.86. (3) 
Combining equations ([2]) and 



rh ^ 18 pc ((t/200 km s"^)^-*^ ^ 13 pc (M./IO^Mq^^^ . (4) 



The M, — a relation extends at least down to M, ^ lO^'^M©, the mass of the Milky Way 



SMBH 



191]. Indirect, but contested, evidence exists for lower mass, intermediate-mass 



black holes in some low-luminosity spheroids, active galaxies, and star clusters 



2Q| 



The connection between SMBHs and nuclear structure is circumstantial but reasonably 
compelling. Observed core radii are ~ a few Vh in the brightest elliptical galaxies, consistent 
with a model in which the cores were created when stars were displaced by a pre-existing 
binary SMBH j21|. At the other extreme in spheroid luminosity, NSCs appear to sometimes 
co-exist with SMBHs, but there are only a handful of galaxies in which the presence of both 
components can unambiguously be established [22|; thus there is no clear evidence that 
SMBHs are associated with an excess of (luminous) mass at the centers of galaxies. 

The nuclear relaxation time 

0.34^3 , , 

tr- = — ^ (5) 



G^Pi^rrii, In A 

measures the time scale over which gravitational encounters redistribute energy between 
stars; InA 12 is the Coulomb logarithm 23|]. Relaxation times greatly exceed 10 Gyr at 
all radii in spheroids more massive than ~ lO^^'M© 2J]; in these "coUisionless" systems, star- 
star and star-dark matter interactions occur too rarely to significantly alter the distribution 
of either component over the lifetime of the galaxy. In fainter spheroids, and particularly 
those containing dense NSCs, central relaxation times can be shorter 25|]; for instance, at 
the Galactic center, tr falls below 10 Gyr inside Vh 26|. In these "collisional" nuclei, the 
distribution of stars around a SMBH is expected to evolve, in a time ~ t^, to the quasi- 
steady-state form 

p^(r) oc r^''/^ (6) 

at r < Th. a Bahcall-Wolf cusp {27]. If multiple mass groups are present, equation ([6]) 
describes the central behavior of the most massive component, while the lowest-mass com- 
ponent (e.g. dark matter particles) obeys 



p(r) oc r ^'^^ 



(7) 



28l | . Equation ([6]) approximately describes the distribution of luminous stars at the Galactic 



center 



26l ] , but no other galaxy containing a SMBH is near enough that a Bahcall-Wolf cusp 



could be resolved even if present. 

Scaling relations between spheroid luminosities and masses or velocity dispersions are 
continuous over many decades in mass, from giant elliptical galaxies down to globular clusters 
(Fig.lH). Only the class of dwarf spheroidal galaxies (dSphs) depart systematically from these 
relations, in the sense of having too large an inferred (dynamical) mass compared with their 
luminosities: these systems appear to be dark-matter dominated even at their centers (§ rVI|) . 



0.5 1 2 5 10 1 2 2 -2 -4 -6 
Sersic n log [kpc] log(Pg) 

FIG. 1: Relations between the parameters that define the mass, size and density of galaxy spheroids 

and A^-body dark matter halos. Because both types of system are well described by the same 

empirical density law, they can be plotted together on the same axes, n = Sersic index, = 

effective (projected half-mass or half-light) radius, pe = mass density at r = R^. The filled symbols 

are dwarf and giant elliptical galaxies; structural parameters for these objects were derived by fitting 

Sersic laws, eq. ([T]), to their observed luminosity profiles. Luminosity densities were converted to 

mass densities using an assumed value for the stellar mass-to-light ratio, neglecting dark matter. 

n 

The open (red) stars are a set of simulated, galaxy-sized dark matter halos [fi]. Mass density 
profiles for these objects were fit with de-projected Sersic laws and total mass was defined as mass 
within the virial radius. The solid lines are semi-empirical fitting relations. (Adapted from 



II. DARK MATTER MODELS 



A. Cusps vs. cores 

Traditionally there have been two approaches to estimating the density of dark matter 
at the centers of galaxies. Unfortunately, they often lead to different conclusions. 

N-body simulations of gravitational clustering follow the growth of dark matter halos as 
they evolve via mergers in an expanding, cold-dark-matter universe. Halo density profiles in 
these simulations are well determined on scales 10^^ < r/vmr < 10°, where the virial radius 
is of ordcr 10^ kpc for a galaxy like the Milky Way; hence inferences about the dark 
matter density on parsec or sub-parsec scales require an extrapolation from the A^-body 
results. A standard parametrization of pdm in these simulated halos is 

PDM(r)=por'(l + 0"' (8) 
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[2!, the NFW profile, where ^ = r/rg and is a scale length of order r^r- In the Milky Way, 
r^r ^ Rq (the radius of the Solar circle) hence eq. ([8]) is essentially a power law at r < Rq 
and the implied dark matter density is 

-1 



PDM(r) ^ IO'MqPc" 



Rcr 



© 



(9) 



1O-2M0PC-3 y V8 kpcy VI pc 

where = Pdm(-R©) and P0 ~ 8 x 10~^MqPc~^ (from the Galactic rotation curve). 

Rotation- curve studies of low-surface-brightness spiral galaxies are generally interpreted 
as implying much lower, central dark matter densities 30|-|34| . While there are caveats to 



';his interpretation - systematic biases in long-slit observations 



35j, non-circular motions 



36|, gas pressure [37|, etc. - these effects do not seem capable of fn 



discrepancies between rotation curve data and expressions like (IE]) 
Pdm('") that is often fit to rotation curve data is 



33 



ly explaining the 
A model for 



PDM[r) 



(10) 



the Burkert profile [30], where C, = r/rc and Tc is the core radius. Inferred core radii are 
~ 10^ — 10^ pc and inferred central densities are 10^^ ^ Pc ^ lO^Mopc"^. 

Since the A^-body halos are not resolved on the scales (~ 10^ pc) where rotation curves are 
measured, the mismatch between theory and observation may be due in part to a poor choice 
of empirical function used to describe the iV-body models. An alternative parametrization 



pDuir) = Poexp - (r/ro)^''" 



(11) 



the Einasto profile 



eq 



39], has recently been shown to describe A^-body haloes even better than 



40l |. The low central density of the Einasto model alleviates some, but not all. 



of the disagreement with rotation curve studies [29 1. 

Remarkably, eq. ( iTTl) has the same functional form as Sersic's law ([I]) that describes the 
projected density profiles of galactic spheroids. In fact, the two descriptions are roughly 
equivalent if nsinasto ~ ""-Sersic + 1 showing that luminous spheroids and simulated dark- 
matter halos are essentially rescaled versions of each other (Fig. [1]), at least over the range 
in radii that is resolvable by the A^-body simulations. 



Several resolutions have been suggested 



measured, central dark matter densities 



41 



br the persistent confiict between predicted and 



43( 1 ■ but none is universally agreed upon. 
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B. Effects of baryonic dissipation 



A^-body simulations of dark matter clustering typically ignore the influence of the baryons 
(stars, gas) even though these components may dominate the gravitational force in the inner 
kiloparsec or so One simple, though idealized, way to account for the effect of the baryons on 
the dark matter is via adiabatic contraction models, which posit that the baryons contracted 
quasi-statically and symmetrically within the pre-existing dark matter halo, pulling in the 
dark matter and increasing its density in the process |44]. When applied to a dark matter 
halo with the density law ([8]), i.e. pdm ~ the result is the more steeply rising p ~ r'^", 

3- 



1.5 



45 



III. DARK MATTER IN COLLISIONLESS NUCLEI 

The dark matter annihilation signal from a region of volume V is proportional to (Pdm)^- 
If the dark matter density rises steeply toward the center of a galaxy, the annihilation flux can 
be dominated by dark matter within the central parsec or so. Neither A^-body simulations, 
nor rotation curve studies, are a reliable guide to pdm on these small scales. In addition, in 
many galaxies, the total gravitational force in the inner parsecs is dominated by the SMBH. 

We consider first "coUisionless" nuclei, in which central relaxation times exceed ~ 10 



Gyr; this is the case in spheroids more massive than ~ 10^^ Mq 2^- In these systems, 
the distribution of stars and dark matter near the galaxy center has probably remained 
essentially unchanged since the era at which the nucleus and the SMBH were created. 

A. Black hole adiabatic growth models 

If the SMBH grew to its final size in the simplest possible way - via spherically-symmetric 



infall of gas - the density of matter around it would increase [49|, |50|, in the same way 
that contracting baryons steepen the dark matter density profile on somewhat larger scales 
(§ IIIBp . In the limit that the growth timescale of the SMBH is long compared with orbital 
periods, this scenario predicts a final density (of stars or dark matter) near the SMBH of 

Pfir) ^ P^{rf){r/rfr\ 7/ = 2 + 1/(4 - 7.) (12) 
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FIG. 2: Possible models for the dark matter distribution at the center of a galaxy like the Milky 
Way. The thin curve shows an Einasto density profile, eq. (fTTI) . and the thick curve labelled "cusp" 
is the same model after "adiabatic compression" by the baryons (stars and gas). Lines labelled 
"spike" show the additional increase in density that would result from spherically-symmetric growth 
of the SMBH. The "annihilation plateau" is the density that satisfies p = m^/ {av)t\ this density 
was computed assuming = 200 GeV, {av) = 10^^^ cm^ s^^ and t = 10^^ yr. Dashed vertical 
lines indicate the SMBH's Schwarzschild radius (left) and gravitational infiuence radius (right). 
Effects of the dynamical processes discussed in ^IVI (scattering of dark matter off stars, loss of 
dark matter into the SMBH) are excluded from this plot; in nuclei (like that of the Milky Way) 
with relaxation times less than ^ 10 Gyr, these processes would generally act to decrease the dark 



matter density 



be. 



(Adapted from 48|].) 



ow what is shown here, particularly in the models with a spike (cf. Fig. I4l5p . 



where pi oc r~"'^ is the pre-existing density profile, and VfK. 0.2r/i. Even for 7^ 0, eq. f ll2p 
predicts 7/ > 2 - a density spike (Fig. |2]). Such a steep dark matter density profile near the 
SMBH would imply very high rates of dark matter annihilation |51 |. 

Stars would respond in the same way as dark matter particles to the growth of a SMBH. 
A p ~ density cusp in the stars is not observed aX r < in any galaxy however, even 
those close enough that a spike could be resolved if present. In the case of the most luminous 
galaxies, this is an expected consequence of core formation by binary SMBHs, as discussed 
in the next section. In low-luminosity spheroids like the bulge of the Milky Way, relaxation 
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times are short enough to convert a stellar spike into a shallower, Bahcall-Wolf cusp in ~ 
one relaxation time. 

It is also possible that spikes never form. Even small (compared with r^) and temporary 
displacements of the SMBH from its central location are sufficient to inhibit the growth of a 
spike or to destroy it after it has formed 52(]. Most models for the growth of SMBHs invoke 
strong departures from spherical symmetry during galaxy mergers in order to remove excess 
angular momentum from the infalling gas 53 1. 



B. Binary black holes and core creation 



Most spheroids are believed to have experienced at least one major merger (defined as a 
merger with mass ratio ~ 3 : 1 or less) since the epoch at which the SMBHs formed 55|. 

n 

If two merging galaxies each contain a SMBH, a massive binary forms pu\ , displacing stars 
and dark matter as the two holes spiral in to the center. The process can be understood as a 
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3: The effect of mergers, including SMBHs, on the central densities of dark matter halos 
54l |. The curves labelled M and S are the density profiles of the larger of the two haloes before 
and after a SMBH was grown adiabatically at its center. The other curves show the final density 
profile of the merged halo, for various choices of the initial halo parameters, as given in the table. 
Ml (M2) is the mass of the large (small) dark matter halo and p is the central halo density before 
growth of the SMBH. Mass and length scalings were based on the Milky Way. 
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kind of dynamical friction, with the "heavy particles" (the SMBHs) transferring their kinetic 
energy to the "light particles" (stars, dark matter). However, most of the energy transfer 
takes place after the two SMBHs have come within each other's spheres of influence, and in 
this regime the interaction with the background is dominated by another mechanism, the 
gravitational slingshot js^]- The massive binary ejects passing stars or dark matter particles 



at high ve 



energy 



58| 



ocity, removing them from the nucleus and simultaneously increasing its binding 



This process stops, or at least slows, when the two SMBHs reach a separation rstaii ~ 
g/(l + g)^r/i, the stalling radius; here q = m2/mi < 1 is the binary mass ratio. At this 
separation, the binary has already removed essentially all material on intersecting orbits 
and the inspiral stops; or, it continues at a much lower rate that is limited by how fast the 
depleted orbits can be repopulated [2J]. The size of the low-density core that is produced 



by inspiral from r ^ to r ^ rgtaii is ~ a few times Vh, quite consistent with the sizes of the 



stellar cores observed in many galaxies 2l|, |59|] . Dark matter cores would presumably be of 



similar size (Fig. [3]), or even larger if multiple mergers occurred 2l|, |60[ or if the pre-binary 
dark matter distribution was characterized by a core as in the Burkert model described 
above. This mechanism can probably not explain the kpc-scale dark matter cores inferred in 
many spiral galaxies, however: the mergers that formed the bulges of these systems would 
have resulted in much smaller, parsec-scale cores. 

C. Gravitational- wave recoil 

If the two SMBHs at the center of a merged galaxy manage to overcome the "final- 
parsec problem" and coalesce, another mechanism comes into play that can affect the central 
density of stars and dark matter. Emission of gravitational waves during the final plunge is 
generically anisotropic, resulting in a transfer of linear momentum to the coalesced SMBH 



6l| . The resultant "kick" can be as large as ~ 4000 km s ^ if the two holes have equal mass 



and optimal spins (i.e. maximal amplitude, oppositely aligned, and parallel to the binary 



orbital plane) 62|. While such extreme kicks are probably rare, even a mass ratio of 0.1 
can result in kicks of ~ 1000 km s~^ if spins are optimal; while if the spins are maximal but 
oriented parallel to the orbital angular momentum, the kick velocity peaks at ~ 600 km s^^ 
for mi = 7712- By comparison, kicks large enough to remove SMBHs from galaxy cores range 
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from ~ 90 km s ^ for spheroid masses of 3 x lO^M© to ~ 750 km s ^ for Mgph = 3 x lO^^M© 
to ~ 1000 km s-i for M^ph = 3 x 10^^ Mq gJ- 

Sudden removal of the SMBH from the gala xy center impulsively reduces the force that 



binds stars and dark matter to the center 



63 



64l |. If Vkick is less than the galaxy central 



escape velocity, still more energy is injected into the core by the kicked SMBH as it passes 
repeatedly through the center before finally coming to rest. Cores enlarged in this way can 
be several times larg er than rh, and indeed a few of the brightest elliptical galaxies have such 
over-sized cores [65| ; dark matter cores are presumably of comparable size in these galaxies. 



IV. DARK MATTER IN COLLISIONAL NUCLEI 



Nuclear relaxation times fall below 10 Gyr in spheroids fainter than ~ 10^^ Lq, roughly the 
luminosity where NSCs first appear 2J]. As discussed above, at least some of these galaxies 
(including the Milky Way) also contain SMBHs. In these collisional nuclei, a Bahcall-Wolf 
cusp in the stars can re-form even if it had been previously destroyed by a binary SMBH 



661]. Both the Milky Way and the nearby dwarf elliptical galaxy M32 exhibit steeply-rising 



stellar density profiles within the influence radii (r < 1 pc) of their SMBHs 



15 



67| 



Dark matter particles in these galaxies are still collisionless: their individual masses are so 
small that gravitational encounters between them are negligible. But even massless particles 
can scatter off of stars, and the associated time scale is equal to within a factor of order 
unity to the star-star relaxation time, eq. (j 
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68 



69|. 



Naively, one would expect the stars to act like a heat source, transferring kinetic energy 
to the dark matter particles and lowering their density. This does occur; but in addition, the 
phase-space density of dark matter particles is driven toward a constant value as a function 
of orbital energy, f{E) ~ /q. A constant phase-space density with respect to E implies a 
configuration-space density that rises as p ~ ^-3/2 ^j-^g potential of a SMBH. The 
term crest, for "coUisionally-regenerated structure" has been coined to describe the result 
of this process [TOj. 

Fig. m illustrates the joint evolution of the stellar and dark- matter densities near a SMBH 
at the center of a galaxy in which the density of both components had previously been 
lowered by a binary SMBH. The stars are seen to attain the p ~ r"'^/'^, Bahcall-Wolf form in 
approximately one relaxation time. Simultaneously, the dark matter particles evolve to the 
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FIG. 4: Joint evolution of the stellar (left) and dark-matter (right) densities around a SMBH due to 



star-star and star-dark matter gravitational encounters [TOj]. Length unit rh is the SMBH influence 
radius (roughly 1 pc at the Galactic center). Density is in units of its initial value at rh- Curves 
show density profiles at times (0, 0.2, 0.4, 0.6, 0.8, 1.0) in units of the initial relaxation time (Eq. [5]) 
at r = r/j. Dashed lines are the "steady-state" solutions, eqs. ([6]) and d?]). 

shallower p ~ r"^/^ profile, increasing their density near the SMBH. The normalization of 
the dark matter density continuously drops as the stars transfer kinetic energy to the dark 
matter. Simulations like these suggest that the presence of a Bahcall-Wolf (coUisional) cusp 
in the stars at the center of a galaxy should always be associated with a shallower^ ~ ^-3/2 



"crest" in the dark matter, regardless of how the nucleus and the SMBH formed [70| . Note, 
however, that this argument can not be used to infer the normalization of the dark matter 
density. 

Conditions for the formation of crests are relaxed somewhat if there is a top-heavy 
spectrum of stellar masses since the dark matter scattering time scales as rh~^ where 
= {ml)/{m^) is the second moment of the stellar mass function 69|]. The stellar cusp 



can also evolve more quickly in this case [71]. However, due to their gradual dissolution, 
dark matter crests might only be present, with significant amplitudes, in galactic nuclei hav- 
ing a fairly narrow range of properties: older than ~one relaxation time but younger than 
many relaxation times. The corresponding range in spheroid luminosities is approximately 
3 X IO^Lq ^ L <3 X 10^ Lq [tqI. In addition, low-luminosity spheroids may not all contain 
massive black holes. In the absence of a SMBH, the stars would undergo core collapse, in a 
time that shorter than < 10^° yr in the densest nuclei |25|, producing a r~^'^^ stellar density 
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profile; the dark matter density would be expected to evolve only slightly in this case 



72| 



V. THE GALACTIC CENTER 



The proximity of the Galactic center makes it a promising target for indirect detection 
studies: predicted fluxes can be more than an order of magnitude higher than for any other 



potential galactic source 
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75| . In addition, observations of stellar velocities in the inner 
parsec of the Milky Way yield a highly precise value for the mass in the SMBH 
as well as (somewhat less precise) estimates of the distributed mass [7g, [TA . In principle, 

nn 

dark matter might be detected by observing its effects on the stellar orbits [78|, |79|, but 
uncertainties about the masses associated with other "dark" components - neutron stars, 
stellar mass black holes, etc. - probably render this approach unfeasible for the forseeable 
future. 

Given a detector with angular acceptance AQ sr, the observed flux of photons produced 
by annihilation of dark matter particles is [t^ 



<^{An,E) ^ 1.9 X 10 
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dN 



(av) 



dE 10-26cm-3s- 



/ ITeV V 



JAnAn cm ^s 



(13) 



where dN/dE is the spectrum of secondary photons per annihilation, is the particle 
mass, (av) is the velocity-averaged self-annihilation cross section, and J^n contains the 
information about the dark matter density: 



J 



AQ 



KAQ- 



An Jtp 



(14) 



Here, dl is an element of length along the line of sight and ip is the angle with respect to the 
Galactic center. The normalizing factor K is typically set to K^^ = (8.5kpc)(0.3 GeV/cm^)^: 
the product of the distance to the Galactic center, and the squared, local value of the dark 
matter density, the latter derived from the measured rotation curve assuming an NFW halo. 
Henceforth we write JAn=io-^ = Jb- AQ = 10^^ sr (~ 10 arc minutes) is the approximate 



angular resolution of atmospheric Cerenkov telescopes like H.E.S.S. 



and of FERMI 



Extrapolation of a halo model like that of eq. (|8]) into the Galactic center region gives 
Js ~ lOf, large enough to produce observable signals for many interesting choices of (erf) 



and m 



74| . On the other hand, the detections by the Whipple and H. 



of 7 rays from the Galactic center with energies up to 10 TeV 



82 



l.S.S. collaborations 



831] would require very 
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FIG. 5: Left: Time evolution of the dark matter density at the center of the Milky Way, including 
the effects of scattering off of stars; self-annihilations; and absorption into the SMBH. The initial 
density profile, shown by the dashed line, represents an adiabatically-compressed NFW profile 
with p 7--2-3_ rpj-^g solid, dotted, and dash-dotted lines show evolved density profiles assuming 
scattering only; annihilation only; and both annihilation and scattering, repectively. The top curve 
in each set is for t = 2.5 Gyr and the bottom curve is for t = 10 Gyr. Right: Evolution of the 
form factor J5 for the models in the left panel. In the absence of self-annihilations (solid line), the 
central density remains very high until ~ 1 relaxation time has elapsed, yielding large J5 values. 
When annihilations are included (dotted, dash-dotted lines), the dark matter density near the 
SMBH drops rapidly and J5 is much smaller at early times. The almost equal values of J5 at 10 
Gyr in the two models that include scattering is partly a coincidence due to the assumed value of 
the annihilation cross section. (Adapted from [s?].) 



large values of {av)J 84j-l86|. motivating the exploration of models in which the central 



dark matter density is enhanced with respect to standard models - for instance, via the 
collisionless "spikes" discussed above. 

In models with such high initial densities, pdm evolves rapidly near the SMBH, due both 
to scattering off of stars and to a number of other processes. 

• Self-annihilations: The same annihilations that produce observable radiation also 
cause the dark matter density to decay. Self-annihilations limit the density to 



Pa 



m 



/{crv)t, with t ^ 10 Gyr the time since formation of the spike SSj. The 



result is a weak, 



-0.5 



density plateau near the SMBH 89|]. Assuming a "max- 



imal" {av) ^ 3 X 10 cm^ s ^, appropriate for a thermal relic, and a "minimal" 
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^ 50 GeV, appropriate to neutralinos in theories with gaugino and sfermion mass 



unification at the GUT scale 



90|, imphes pann ~ lO^M© pc^^ at t = 10 Gyr. 



Capture of dark matter within stars. Another potential loss term for the dark matter 
is capture within stars, due to scattering off nuclei followed by annihilation in stellar 
cores. However this process is not likely to be important unless the cross section for 
WIMP-on-proton scattering is very large. 

Capture of dark matter within the SMBH. Any dark matter particles on orbits that 
intersect the SMBH are lost in a single orbital period. Subsequently, scattering of 
dark matter particles by stars drives a continuous flux of dark matter into the SMBH 



9l|. Changes in orbital angular momentum dominate the flux; in a time ~ t^, most 
of the dark matter within rh will have been lost, although the net change in the dark 
matter density profile will be more modest than this suggests since more particles are 



continuously being scattered onto depleted orbits 



69] 



A strict inner cut-off to the dark matter density is set by the SMBH's event horizon, vsch = 
2GM,/c^ ^ 3 X 10^^ pc although for reasonable values of and (erf), the density is limited 
by self-annihilations well outside of rsch- 

These various effects can be modelled in a time- dependent way via the orbit-averaged 



Fokker-Planck (FP) equation 69|]. In its simplest, energy- dependent form, the FP equation 
can be written 

i = a# - SiEy^E) - mME). (15) 

Here f{E) is the phase-space mass density of dark matter, E = —v'^/2 + 0(r) is the energy 
per unit mass of a dark-matter particle, p is a phase-space volume element, and 0(r) is the 
gravitational potential generated by the stars and the SMBH. Fe is the flux of particles 
in energy space and depends on the stellar density profile and the stellar masses. The two 
loss terms, i^cou and i^ic, represent decay of dark matter due to self-annihilations and capture 
within stars; and loss of particles into the SMBH, respectively. A similar equation describes 



the evolution of the stellar distribution [70|. 



Self-annihilations add a second time scale to the problem that depends on the particle 
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physics parameters, the annihilation time 

lann = PDM 

-nsrvrT ^( ^™ VV ^^""^ V' (m) 

• ^ VlOO GeV/ VIOIOM0PC-3; Vl0-26cm3s-V ' ^ ^ 

Even assuming the "maximal" annihilation model defined above, tann becomes comparable 
to tr only for dark matter densities greater than ~ IO^Mq pc~^: corresponding to r < 10^^ 
pc in models with a "spike." Fig. [5] illustrates this. At early times, annihilations dominate 
the changes in J, but after ~ 1 Gyr heating of dark matter particles by stars tends to 
repopulate orbits near the SMBH, tending toward a pdm ~ r~'^/^ "crest" density profile at 
r < Th- 

The predicted spectrum of annihilation products depends separately on (shape) and 
{av) J (amplitude) 9^ ; while for a given initial dark matter model, the final distribution of 
mass in the evolutionary models depends on tr/^ann oc {av)lm^. Assuming a dark matter 



origin for the TeVgamma rays observed by H.E.S.S. |83| . the spectrum implies lOTeV < 
niy^ < 20TeV and {av)J^ 10~^^ cm^ s~^ 86}. The latter value requires either J5 > 10^, or 
a substantial enhancement in the dark matter relic abundance compared with expectations 
for thermal freeze-out. Fig. [5] suggests that a dark matter "crest" can yield sufficiently high 
values of J if the initial density profile is sufficiently steep. 

VI. DWARF SPHEROIDAL GALAXIES 

As noted above, dwarf spheroidal (dSph) galaxies depart systematically from the scaling 



relations obeyed by other "hot" stellar systems between size, mass and luminosity 93 |. 
dSphs have total luminosities and internal velocity dispersions comparable with those of 
globular star clusters, but are much larger, implying very high ratios of (dynamical) mass to 



(stellar) light, roughly 10 — 100 times the Solar value [9J]. Since the stars in these systems 
contribute a negligible fraction of the total mass, dSphs are ideal test-beds for theories of 
dark matter: in prin ciple, Pdm('") can be mapped directly given sufficiently large samples of 



radial velocities 



95|, without the need to correct for baryon "contamination.'' 



Modelling of this sort has been carried out now for roughly a dozen dSphs [e.g. 



96- 



98] . In practice, the kinematical data are not copious enough for a fully nonparametric 



approach and the inferred densities are still somewhat model-dependent. Interestingly, the 
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kinematical data seem rarely if ever to demand cusped dark matter density profiles like 
;hat of e g . ([8]). Halo models with low-density cores, e.g. eq. flTOj) . are sometimes preferred 



99 



100|. Assuming that /3dm(^) follows an NFW profile, inferred mean densities within 
10 pc are ~ IOMq pc^^; compared with ~ O.IMq pc~^ if the inner den sity profile is flat. 



Inferred total masses depend less strongly on the assumed profile shape 101| . Thus, while 
the existence of large amounts of dark matter is clearly established in the dSph galaxies, 
as expected based on the A^-body models, the detailed distribution of mass within these 
systems appears to be at odds with the A^-body predictions. dSph galaxies are similar to 
low-surface-brightness spiral galaxies in this regard. 

The number of k nown dSph satellites of the Local Group has roughly doubled during 



the last decade 



102| and will probably continue to rise; at last count the Milky Way halo 



contained at least 18 dSphs. Their proximity, c omb i ned with their large masses, make them 
good candidates for indirect detection studies 



103 



104( 1 . althou gh th e predicted fluxes are 



interestingly large only if the dark matter distribution is cusped |105| . 

Angular sizes are small enough that a large fraction of their dark matter could be imaged 
in a single pointing of a telescope like FERMI. The predicted annihilation flux then scales 
simply as ~ p^rf where ps is the dark matter density at the scale radius r^. If the dark 
matter is clumped, fluxes could be boosted by up to two orders of magnitude |106| . In 
addition, the low baryonic densities imply that dSphs should be relatively free of intrinsic 
gamma rays from other astrophysical sources, making the interpretation of a signal much 
more straightforward than in the case of the Galactic center. A particularly attra ctive object 
is the recently-discovered object Willman 1 with a luminosity of only ~ 10^ Lq 107l |. At a 
distance of 38 kpc, this object is so close that it could be marginally resolved by FERMI, 



in principle allowing a determination of the dark matter distribution 



106| 



I think A. Graham and E. Vasiliev for supplying modified versions of figures from their 
published work. 
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